INVERSE PROBLEMS AND INDEX FORMULAE FOR 
DIRAC OPERATORS 



YAROSLAV KURYLEV AND MATTI LASSAS 

Abstract: We consider a Dirac-type operator Dp on a vector bundle V 
over a compact Riemannian manifold (M, g) with a nonempty bound- 
ary. The operator Dp is specified by a boundary condition P{u\qm) — 
where P is a projector which may be a non-local, i.e. a pseudodiffer- 
ential operator. We assume the existence of a chirality operator which 
decomposes L 2 (M, V) into two orthogonal subspaces X + © X_. Under 
certain conditions, the operator Dp restricted to X + and X_ defines 
a pair of Fredholm operators which maps X + — > X_ and X_ — > X + 
correspondingly, giving rise to a superstructure on V . In this pa- 
per we consider the questions of determining the index of Dp and 
the reconstruction of (M,g), V and Dp from the boundary data on 
dM. The data used is either the Cauchy data, i.e. the restrictions 
to dM x R + of the solutions to the hyperbolic Dirac equation, or the 
boundary spectral data, i.e. the set of the eigenvalues and the boundary 
values of the eigenf unctions of Dp. We obtain formulae for the index 
and prove uniqueness results for the inverse boundary value problems. 
We apply the obtained results to the classical Dirac-type operator in 
M x C 4 , M C M 3 . 

AMS classification: 35J25, 58J45. 
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1. Introduction 

Recent study of inverse problems has successfully shown that exter- 
nal measurements can be used to uniquely determine coefficients of vari- 
ous partial differential equations modeling macroscopic and microscopic 
phenomena. Examples of these are the paradigm problem, the inverse 
problem for the conductivity equation encountered in the impedance 
tomography, the inverse problem for the wave equation with an inhomo- 
geneous wave speed, the inverse scattering problem for the Schrodinger 
operator (for different approaches see e.g. pi [HI El EH1 EH HH SHI 152] . 
In these problems the structure of the space is a priori known before 
measurements, being usually a known domain of an Euclidean space. 
Recently, many inverse problems have been generalized to the cases 
where the underlying space is not a priori known but is assumed to be 
some unknown smooth Riemannian manifold. This was followed by the 
change of the attitude from proving the uniqueness in the considered 
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inverse problem to describing the corresponding groups of transforma- 
tions of the unknown object which preserve the measured data, e.g. 
isometries of Riemannian manifolds, gauge equivalences of differential 
operators, etc. From this point of view any anisotropic inverse problem, 
even in a given domain Del", falls into this second category having a 
non-trivial group of transformations which consists of some boundary 
preserving diffeomorphisms of the domain, see e.g. [U 13 E21 ESI EH1 HH] 
and, for a detailed description of this approach, [2~H] . 

A natural task is to generalize this approach further to inverse prob- 
lems on general vector bundles as, in the modern physics, physical phe- 
nomena are often modeled by equations on bundles. Indeed, the vector 
bundles are encountered from the early development of the relativis- 
tic quantum mechanics to the modern quantum field theory. Thus, 
in this paper our aim is to bring the mathematical study of inverse 
problems closer to the modern physics and to study inverse problems 
for one of its basic equations, the Dirac equation on a general vector 
bundle. Assuming the existence of a chirality operator, we consider a 
Dirac equation on an unknown bundle and reconstruct the underlying 
manifold, the bundle structure on it and also the coefficients of the 
equation up to a natural group of transformations, namely, the isome- 
tries of the manifold and bundlemorphisms of the bundle. To this end 
we further advance the techniques developed for the inverse problems 
for Maxwell's equation [Ml 1371 . 

For previous results in inverse problems for the Dirac equation in the 
1-dimensional case, see e.g. [HJ EHJ El] and in the 3-dimensional case, 
see [2U EEl E7J EE] • In the 3-dimensional case all these papers deal with 
perturbations of the canonical Dirac operator D in M 3 or Q C M 3 , 

3 

(1) D u = i 2_, a kdkU + mc 2 aoU, u G L 2 (R 3 ; C 4 ), 

fc=i 

where a u , v = 0, . . . , 3, are the standard Dirac matrices. Thus the 
structure of the underlying space and Dirac bundle was known a priori. 

1.1. Formulation of the problem. In this paper we consider two 
types of inverse boundary value problems for the Dirac-type operator, 
D, acting on sections of a complex vector bundle, V, over a compact, 
connected Riemannian n-manifold (M, g) with a non-empty boundary, 
dM. Roughly speaking (for a rigorous formulation of the problems see 
Subsection 2.4 below), the first inverse problem deals with waves, i.e. 
the solutions to the hyperbolic problem, 

(2) (idt + D)u(x, t) = in MxM+, u\ t=0 = 0. 

Then the data used is the set of the Cauchy data of these waves re- 
stricted to dM x M + , i.e. the set Cq(D) = {u\qmxr + '■ u satisfies (J2J)}- 
The corresponding inverse problem, called the dynamic inverse prob- 
lem, is that of the unique determination of the manifold M, the bundle 
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V, and the operator D when we are given the boundary dM, the restric- 
tion W of the bundle to dM, W = V\qm, and the dynamic boundary 
data C (D). 

The second inverse problem or, more rigorously, a family of inverse 
problems, deals with the self-adjoint Dirac-type operator, 

(3) Dpu = Du; V(D P ) = {ue H\M, V) : P(u\ 9M ) = 0}. 

Here P is an orthoprojector in L 2 [dM, W), given by a zero-order classi- 
cal pseudo-differential operator on dM which, in general, is not a local 
operator. The corresponding data are the boundary spectral data, 

(Afc, (f>k\dM), k — 1,2, ... , 

where <pk are the eigenvalues and normalized eigenfunctions of Dp. 

Throughout the paper, we assume the existence of a chirality oper- 
ator F : V — > V, acting fiber- wise, i.e., F : 7T -1 (a;) — > 7r _1 (x), with 
7r : V — > M being the projector onto the base M. The chirality opera- 
tor satisfies the anti-commutation relation 

(4) DoF + FoD = 

and makes it possible, under some conditions, to split Dp into a pair 

D -( D ? 



of operators Df 



We will derive formulae to express Ind (Dp) in terms of the dynamical 
boundary data. 

The main tool in our study of the inverse problems for the Dirac-type 
equations is the boundary control method, see |H], or, more precisely, 
its geometric version, see (2H] for a detailed exposition. The method 
is based on the finite velocity of the wave propagation and control 
properties of the waves generated by boundary sources. 

The plan of the paper is as follows: In Section 2 we give definitions of 
the Dirac bundles and Dirac-type operators, introduce a decomposition 
of L 2 (M, V) into two channels X + © X_ and formulate rigorously the 
main results of the paper. In Section 3 we describe properties of the 
hyperbolic Dirac-type equation and provide a unique continuation 
result. In Section 4 we derive a formula for the energy of a wave in each 
channel in terms of the boundary data and prove some global and local 
boundary controllability results for the hyperbolic Dirac-type equation. 
Generalized boundary sources and an extension of the time-derivative 
d t are introduced in Section 5. The relation of d t and Dp is studied 
in Section 6 where we obtain index formulae for Dp in terms of the 
boundary data. In Section 7 we prove that the dynamic inverse data 
determine the manifold M and its metric g. In Section 8 we introduce 
focusing sources which generate waves that, at given time, are localized 
at a single point x G M. These waves with different polarizations are 
used in the reconstruction of the bundle V and operator D P from the 
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dynamic inverse data. In Section 9 we study the high-energy asymp- 
totics of the solutions to prove the uniqueness for the inverse boundary 
spectral problem. In Section 10 we apply the obtained results to the 
Dirac-type equation ((T]) in M 3 leading to an analog of the Aharonov- 
Bohm effect. In Appendix we prove, using the construction of Briining 
and Lesch, the self-adjointness of the Dirac-type operator with the local 
boundary condition. 

2. Definitions and main results. 

2.1. Definitions. Here we introduce some basic definitions, results 
and examples about Dirac bundles and operators on them following 
the terminology of [T3 | ITKl IKK] . 

Let (M, g) be a compact connected C°°-smooth Riemannian n-mani- 
fold with a non-empty boundary. Let V be a smooth complex vector 
bundle over M and denote the projection onto the base manifold by 
7T : V — > M. Each fiber 7r _1 (x) is a complex (i-dimensional vector space 
with a Hermitian inner product (■,•) = (■, -) x . We define \cj)\ x = (4>,(p) x 
and denote the smooth sections of V by C°°(M, V). Endomorphisms 
End (V) are fiber-preserving smooth maps, L : V — > V, vrL(0) = ncf) 
that are linear in each fiber. By Clif(M) we denote the Clifford bundle 
over M. This means that, at each x G M, the fiber Clif x (M) is an 
algebra generated by the vectors in T X M with a product • satisfying 
the relation 

(5) v-w + w-v = -2g(v,w) for v, w G T X M C Clif x (M). 

We assume that there is a fiber- wise map 7 : Clif(M) — > End (V) that 
provides a Clifford module structure for V, i.e. 7 gives an action, for 
any x G M, of the algebra Clif x (M) on fibers 7r _1 (x) of V and satisfies 

(6) ( 7 (»A, n) x + (A, -y(y)fj.) x = 0, v G T X .(M), A, /i G n-\x). 

Let V be the Levi-Civita connection in (M,g). It defines a connection 
on Clif (M) that satisfies 

Vx(d • w) = (Vxv)-w + v ■ (Vi«;), 

where X G C°°(M,TM) and v,w G C°°(M, Clif (M)). We denote by 
V also a connection on V that is compatible with (• , • ) and 7, i.e., 

X(0,^) = (Vx0,^) + (0,V^}, 
Vx(7(«>)0) = 7(Vx^)0 + 7(^)Vx0, 

where X G C°°(M, TM), w G C°°(M, Clif (M)), and 0, ip G C°°(M, 1/). 
When these conditions are satisfied, we say that (V, ( •, -),7,V) is a 
Dirac bundle over (M,g), for brevity, a Dirac bundle V. 
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The unperturbed Dirac operator, D on the Dirac bundle V, which 
is sometimes called the Dirac operator associated with a Clifford con- 
nection, is locally given by 

n 

(7) D Q u\u = ^2^(e j )V ej u\ Ul 

3=1 

where (ei, . . . , e n ) is an orthonormal frame in an open set U C M. 

As noted before, we assume that there is a chirality operator F : 
V — > V acting fiber-wise and satisfying 

(8) (FA, fi) x = (A, Ffi)*, A, fi G 7r _1 (x), 

j(v) o F + F o 1 {v) =0, f G T X M, F 2 = I, 

VF = 0, i.e. Vx(Fu) = F{V x u) 

where X G C°°(M,TM) and u G C7°°(M,V). When n = dim(M) is 
even, a chirality operator always exists, namely, 

(9) F^u = (v /Z T)" /2 7(ei)7(e 2 ) • • • l{e n )4>\u. 

A chirality operator exists also when we deal with a space-type hyper- 
surface on a Lorentzian manifold, e.g. [23j. 

Using a chirality operator, we define the fiber-wise orthogonal pro- 
jectors onto the +1 and —1 eigenspaces of F, 

(10) n+ = i(/ + n n_ = i(/-n 

We denote by Q G End (V) a self-adjoint potential, i.e., 

(11) (Q\n) x = (X,Qff) x , X,fi G 7r _1 (a;), 
which respects the chirality structure, i.e. 

(12) FoQ + QoF = 0. 

We call the minimal Dirac-type operator, D min , the operator 

(13) D min u = D u + Qu, uE V(D min ) = C °°(M, V), 

where Cq°(M, V) is the space of the smooth compactly supported sec- 
tion of V . 

Example 1. As an example, we consider the so-called form-Dirac 
operator, see e.g. [El EE]- Let (M,g) be an oriented Riemannian 
manifold of dimension n and Q?M be the space of the complex- valued 
differential j-forms on M. Let QM = Q°M © • • • © Q n M be the Grass- 
mannian bundle on M and * the Hodge operator of Q(M) related to 
the metric g. For A G Cl 3 x M and \i G fi^M, let (A, \i) x = if j ^ k and 
(A, ji) x = *(A A *jx) for j = k. 

Denote by I : T X M — > T^M the identification I(a 3 -Jj^j) = gjko)dx k 
and by t„ : f2£M — ► f^ _:L M the inner product with a vector t> G T X M, 
e.g. 
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for any w G Q J X M, v,Vi,...,, Vj-i G T X M. Then the map 7 : T X M — > 
End(ft x M), 

7 (u)A = -/(f) A A + t B A 

extends to a homomorphism 7 : Clif^M) — > End(Q x M) . Let V be 
the Levi-Civita connection on QM. Then (QM, ( -, 7,7, V) is a Dirac 
bundle. The unperturbed Dirac operator on this bundle is 

D -MM -> OM, D = d + 5 : f^M -»■ fi J+1 M © fi^M, 

where d is the exterior differential and 

6 = * d* : fi J M fi^M 

is the codifferential. A chirality operator, F for the form-Dirac operator 
may be defined simply as 

(14) FA = (-1) J A, A G Q j M, 

splitting the Grassmannian bundle into the forms of the even and odd 
orders, QM = Q e M © Q°M. 

Example 2. Consider Maxwell's equations in M x R, M C R 3 , 

cmlE(x,t) = —B t (x,t), D(x,t) = e(x)E(x,t), 

cwlH(x,t)= D t (x,t), B(x,t) = fi(x)H(x,t) 

where e(x) and n(x) are positive definite matrix valued functions. The 
velocity of the wave propagation is independent of the wave polar- 
ization if and only if fi(x) = a(x) 2 jj,(x) with some scalar function a. 
In this case the travel time is determined by the Riemannian metric 
gij = a~ 2 det(e)~ 1 5 lk e 3 k . As shown in [HEIEE], Maxwell's system may be 
extended to the hyperbolic equation 

(15) id t uo + i(d - 5 a )uo = 0, lu(x, t) = (uj°, u 1 , u 2 , u 3 ) C QM, 

where 5 a uP = a,8{aT 1 uji) and 5 is the codifferential with respect to the 
travel time metric g. Indeed, a solution of Maxwell's system gives rize 
to a solution of the Dirac-type equation (fl~5|) with u° = 0, u l = Ejdx^ , 
J 1 = B 3 dx l A dx 2 - B 2 dx 1 A dx 3 + B x dx 2 A dx 3 , and u 3 = 0. Defining 
the inner product 

(cu j , X k ) x = aix)- 1 * (u j A *X k )5 jk , for u j G Q j x M, \ k G Q k x M, 

the operator i(d — S a ) is formally self-adjoint in QM. This oper- 
ator defines a Dirac-type bundle (QM, (•,■)%, 7, V) where j(v)X = 
—i(I(v) A A + l v X) and V is the Levi-Civita connection in (M,g). The 
chirality operator jHj) is a chirality operator on this bundle, too. 

2.2. Self-adjoint operators. We return now to the discussion of the 
boundary conditions used to extend the minimal Dirac operator (|13|) 
to a self-adjoint operator in L 2 (M, V) with the inner product 

««»= / (<f>(x),i>(x)) x dV g (x), 

J M 
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where dV g is the Riemannian volume on (M, g). Self-adjoint extensions 
of D min defined by non-local boundary conditions, go back to the cen- 
tennial work by Atiyah-Patodi-Singer p]. In this paper, we will work 
with the self-adjoint extensions of D min extensively studied by Briining 
and Lesch jTHJ , see also [16] f22j. 

In the following we denote by H S (M, V) the Sobolev spaces of sec- 
tions of V with components in the Sobolev spaces H S (U) in local 
smooth trivializations $t/ : 7r _1 (C/) -> UxC d , U C M and by Hq(M, V) 
the closure of C °°(M, V) in H S (M, V). 

Using the natural embedding j : dM — > M , we introduce the induced 
bundle W = V\ dM = fV on dM. 

Next we consider the possible extensions, in L 2 (M, V), of the mini- 
mal Dirac-type operator (fl3|) . We use D to denote the maximal exten- 
sion of the operator -D m in to H l (M, V), 

(16) Du = (A) + Q)u, V(D) = H\M,V). 

It is shown in [12] that there is a wide class of self-adjoint extensions 
of -D m i n defined by boundary conditions of form (JBJ) , 

(17) Dpu = Du, V{D P ) = {ue H\M, V) : P{u\ 9M ) = 0}, 

so that Dp C D. Here P is a zero-order classical pseudo-differential 
operator defining an orthoprojector on L 2 (dM,W). It should satisfy 
the following conditions: 

(18) P 1 (N)= 1 (N)(I-P), 
and 

(19) {P, Paps} is a Fredholm pair, 

that is, Paps '■ R an (P) ~^ R an (Paps) is a Fredholm operator. Here 
N stands for the unit interior normal field of dM and the operator 
Paps is associated with the Atiyah-Patodi-Singer boundary condition. 
Namely, let A(0), sometimes called the hyper surface Dirac operator, 
(see e.g. Appendix or [22].). be a self-adjoint operator in L 2 (dM, W), 

n—l j 

(20) A(U) = - 1 (N)Y,l(e a )V a + r ^-H(x'), V a = V £q . 

o=l 

Here the vector fields e a , a = 1, . . . ,n — 1, form a local orthonormal 
frame on dM, while H{x') is the mean curvature of dM at the point 
a;', where x' — (x\, . . . , x n _i) are local coordinates on dM. Then Paps 
is a spectral projector of A(0) such that 

(21) -P(o,oo) C Paps C -P[o,oo), 

(22) P AP5 = - 7 (iV)(J - PapsHN), 

where we denote by P( a ,b) the spectral projector of A(0) associated with 
an interval (a, b). 



8 



YAROSLAV KURYLEV AND MATTI LASSAS 



Theorem 2.1 (Bruning-Lesch). Let P be a zero-order classical pseudo- 
differential operator defining an orthoprojector which satisfies condi- 
tions $18\) and (1 1 .9)). Let Dp be the Dirac-type operator (Ii7j). Then Dp 
is self-adjoint, the spectrum of Dp is discrete, and all eigenspaces are 
finite dimensional. Moreover, Dp is regular, that is, Dpu G H S (M, V) 
implies that u G H S+1 (M, V) for any s > 0. 

The proof of Theorem 12 .11 is given in |13[ Thm. 1.5] for more general 
Dirac-type operators than those considered in this paper and in [Tfijl^ 
for more restricted cases. 

In the sequel we will denote the eigenvalues of Dp, numerated ac- 
cording to the their multiplicity, by Af , |Af | < |Af +1 |, and the corre- 
sponding orthonormal eigenfunctions by 0{\ When there is no danger 
of confusion we will skip the index P in these notations. 

Using projections 11+ and see (JTU1) we define the bundles V + = 
Ii + V and V- = H-V over M with the projections to M denoted by 
7r±. Note that, for any x G M, tt+ x (x) and 7rZ 1 (x) are orthogonal 
so that V = V + © V- and L 2 {M,V) = L 2 {M,U + V) © L 2 {M, IL_V). 
With a slight abuse of notation we use n ± for the orthoprojectors 
in L 2 (M,V) onto L 2 (M,U±V). The maximal Dirac-type operator D 
satisfies DF + FD = so that D can be decomposed as 

(23) D :H\M, U+V) © H\M, U^V)^L 2 (M, T1_V) © L 2 (M, U+V), 
When 

(24) U + V(D P ) C V{D P ) or, equivalently, n_ V{D P ) C V(D P ), 

decomposition (|2H|l gives rize to a decomposition of Dp. Namely, for 
D% = D P : U ± V{D P ) -> L 2 {M, U T V), 

(25) Dp : U + V(D P ) © UJD{D P ) -> L 2 (M, n + U) © L 2 (M, n_F), 

If (J21j) is valid and thus the decomposition (f23| is possible, we say 
that the bundle V = V + © V_ on M has a superstructure. It then 
follows from [T3| that the resulting operators -Dp are Fredholm with, 
in general, non-trivial indices. For important relations of the index and 
the geometry of the bundle, see e.g. [T| l9| fl8|. fT9]. 

2.3. Local boundary condition. As noted in Introduction, the 
boundary projection P in (fTTj) is usually non-local. However, for a 
Dirac-type operator with chirality, it is possible to introduce a local 
boundary condition which makes it self-adjoint. To this end, we con- 
sider a fiberwise operator T on W = V\qm, 

(26) r : W -> W, TX = Fj(N) A, AG vt^x), x G dM. 
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It is clear from properties © and (jHJ that 

(27) r 2 = j, roF + For = o, 

(T\,fi) x = (\,Tn) x , A,/i G 7r _1 (x). 

Let 

(28) P r =~(/ + r), 

be an orthoprojector acting fiber- wise on 7r _1 (x), x G <9M. Then Pr 
satisfies conditions (|19|) and (fl8|) and, therefore, defines a self-adjoint 
Dirac-type operator D r := D Pr . We note that this operator will play 
a crucial role in the sequel. Because of this, we give in the Appendix a 
sketch of the proof of its self-adjointness. 

Example 1, continued. Clearly, Pr defines a local boundary condi- 
tion for the form-Dirac operator d + 5 determined in Example 1. In 
addition, there are other local boundary condition closely related to 
the de Rham complexes over (M,g). Recall that j : dM — M is a 
natural embedding. The relative boundary condition is defined by 

(29) P r ( u \ 9M ) :=f(u) = 0, 
and the absolute boundary condition by 

(30) P a (u\d M ) :=3*{ln{u)) = 0. 

In the future we denote the corresponding self-adjoint form-Dirac op- 
erators by (d + 6) r and (d + 5) a . 

Example 3. Another important example of a Dirac operator on the 
Grassmannian bundle QM, for n = 2k, is given by 

D x u = (d+5)u, u G Dpi) = {u G H\M,V) : P A ps(u\ 9m ) = 0}, 

where Paps is the Atiyah-Patodi-Singer boundary condition defined in 
subsection 2.2, see e.g. [UHHUS]- This operator is called the signature 
operator. In this case the chirality operator is given by (jOJ. 

2.4. Boundary data for inverse problems. 

Definition 2.2. Consider a Dirac bundle (V, ( •, •), 7, V) with chirality 
F. The induced bundle structure on DM is the collection 

(31) {dM, g\ 9M , W,(-,- ) x \ xed M, l(N)\ 0M , F\ dM }. 

Here dM is considered as a Riemannian manifold with the differen- 
tiable structure induced by the embedding j : dM — > M and the 
metric g\dM = 3*9- The bundle W = V\qm is the induced bundle on 
dM, (■ , • ) x the Hermitian structure to W, j(N)\dM and F\qm are the 
restrictions of the Clifford action j(N) and of the chirality operator F 
on W. 
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Definition 2.3. Let D P be a self-adjoint Dirac-type operator of form 
(fTT|) on a Dirac bundle (V, ( •, -),7, V) with chirality F. The set 

(32) {(A fc , 4>k\dM)}kLi, 

where and 0& are the eigenvalues and the orthonormal eigenfunctions 
of Dp is called the boundary spectral data of D P . 

To define another type of boundary data, observe that any Dirac- 
type operator (fTTjl is associated with a hyperbolic initial-boundary 
value problem 

(33) (id t + D)u(x, t) = in M x E+, 

(34) P(u\ aMxK+ ) = fe PC°°(dM x R + , WO), u\ t=0 = 0, 

o 

where C°°(dMxR +1 X) is the class of the smooth functions on dMxR + 
(with values in X) which are equal to near t — 0. We denote by 
u = w the solution of (f33l - (f34l) . Problem (f33l) - (f34l) defines the response 
operator, Ap, 

(35) hpf = U f \ dM xR + - 

When P = Pp, we denote A r = A Pr . Other data for the inverse 
boundary problem could be the response operator Ap. However, we 
prefer to work with data which is independent of P. 

Definition 2.4. Let D be the maximal Dirac-type operator (|T6|) . Then 
the Cauchy data set of D is the set 

(36) C (D) = {u\ dMx R + : u E C°°(M x R+, V) satisfies ©}. 

Later we prove the following equivalence of different type of data. 

Lemma 2.5. Assume that the induced bundle structure ^Sl\) on dM 
is given. Then the Cauchy data set Cq(D) determines the map Ap and 
vice versa. Moreover, the set Cq(D) and a projector P satisfying $18\) 
and JZi|) determine the map Ap and vice versa. 

2.5. Main results. First we formulate our main index formula. 

Theorem 2.6. Let V be a Dirac bundle over (M,g) with a chirality 
operator F. Let Dp be a self-adjoint Dirac-type operator of the form 

Assume that we are given the induced bundle structure on dM 
and the response operator Ap. Then these data determine, in a con- 
structive way, a Hilbert space T , a pair of projectors B + and B- in T 
with B + @B_ = I , and an unbounded self-adjoint operator d t : T — > T 
with T>(d t ) C T such that 

1. The condition fH^ ) for Dp is valid if and only if B + V(d t ) C T>(d t ). 

2. If decomposition \2J$ is valid then the operator 

(37) Of : B + V(d t ) -> B_T, d+u = d t u 
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is a Fredholm operator and 

(38) Ind(dt) = Ind(D%). 

In Theorem 12.61 the Hilbert space T is the completion of the func- 
tions PC °°(<9M x [0,T],W)), where T > is sufficiently large. This 
completion is taken with respect to a seminorm || • that is explicitly 
determined by the Cauchy data set C (D) (see (|50|) . (|5T|) . (|73|) below). 
The operator d t in T is an extension of the time derivative 

d t f = ^, forfePCSr{dMx[0,T\,W)). 

Next we formulate our main results on the reconstruction of the 
manifold and the Dirac bundle. To this end, let V and V be vector 
bundles over diffeomorphic manifolds M and M. Let i : M — > M be a 
diffeomorphism and it : V — > M and n : V — > M be projections on the 
base spaces. A map L : V — > V is called a bundlemorphism compatible 
with £ if £(vr0) = n(L<f)), <p G V and 

L : 7r _1 (x) — > tt~ (£(x)) is a linear isomorphism for anyx G M. 

Dealing with Riemannian manifolds and Dirac bundles on them we 
need to extend the notation of a bundlemorphism. 

Definition 2.7. Let (V, (■,■), 7, V) and (V, ( ■, ■ V) be Dirac bun- 
dles over isometric Riemannian manifolds (M, g) and (M , g) with chi- 
rality operators F and F, and let D and .D be Dirac-type operators of 
form (fl6jl on them. A bundlemorphism L compatible with an isometry 
£ : M — > M is called a Dirac bundlemorphism if 

(39) (L(p,Lip) e{x ) = 0,^ € 7r _1 (^), 
7(^)) =L 7 ( W )L- 1 , d G Ta-M, 

F = LFL' 1 , D = LDL-\ 

where d£ is the differential map of £. 

We turn now to some induced structures on the boundary. 

Definition 2.8. Let (V, ( -, -),7,V) and (V, ( ■ V) be Dirac bun- 
dles over manifolds M and M with = V\qm and = V\qm- Let 
<9M and 9M be isometric with an isometry k : DM — > 9M. We say 
that a bundlemorphism F" : — > is an isomorphism between the 
induced bundle structures of <9M and <9M compatible with K if 

Kip) K ( x ) = (0, V),, F„ (sb) = KF x K-\ 

7(tZ/c(v)) = ^(v)^" 1 , 7(AT(«(a;))) = ^(iV^iT 1 . 

Here x G <9M, G 9M, 0,^ G vr" 1 ^) and u G T x (dM). At last, 
iV(x) and N(x) are the unit interior normal vectors to dM at x and to 
<9M at x. 
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Let now u G C°°(M x R + , V) be a solution of the hyperbolic Dirac 
equation (|33|) and L : V — > V be a Dirac bundlemorphism. Then 
u = Lu is a solution of the hyperbolic Dirac equation ijHHJ) with D 
replaced with D. Therefore, the Cauchy data sets Cq(D), Cq{D) satisfy 

K{C {D)) = C {D), 

with K = L\w- This shows that the following uniqueness result is 
optimal. 

Theorem 2.9. Let D and D be Dirac-type operators of the form ITflJ) 
on the Dirac bundles with chirality (V, F) and (V, F). Assume that 
there is an isometry k : dM — > dM and an isomorphism between the 
induced bundle structures of dM and dM , K : W — > W, which is 
compatible with k. Assume also that the Cauchy data sets satisfy 

K{C {D)) = C (D). 

Then there is an isometry £ : (M,g) — > (M,g) with £\om — K and 
a Dirac bundlemorphism L : V — > V with L\w = K. In particular, 
D = LDL L . 

Note that this theorem states that the Dirac operators are equal up 
to a bundlemorphism, rather than the connections V and V. 

For the inverse boundary spectral problem we obtain a similar result. 

Theorem 2.10. Let Dp and Dp be self-adjoint Dirac-type operators 
of form ( TTTj ) on the Dirac bundles with chirality (V, F) and (V, F). 
Assume that there is an isometry k : dM — > dM, an isomorphism 
K : W — > W compatible with k, and that the eigenvalues and the 
normalized eigenf unctions of Dp and Dp satisfy 

\j = \j, (f)j o k\ 9m = K{<j)j\ aM ), j G Z + . 

Then there is an isometry £ : (M,g) — > (M, If) with £\om = k and a 
Dirac bundlemorphism L : V — > V with L\w = K such that 

Dp = LD P L-\ P = KPL-\ 

Remark 2.11. To consider above results in more practical terms, as- 
sume that we are given the boundary dM, the bundle W on it, and 
the Cauchy data set C (D) or the boundary spectral data. The proofs 
of the above theorems are constructive, providing a method to recover 
the manifold M (up to an isometry), the bundle V with its Hermitian, 
(•, •) and Clifford module, 7 structures, and the operators D and F 
(up to a Dirac bundlemorphism). 
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3. Hyperbolic Dirac equation and unique continuation 

3.1. Basic properties. We start with a general boundary condition 
P satisfying (|18| and (|19|) . Denote by u^(t) = w(x,t) the solution of 
the initial boundary value problem (j33| - (|34j) . 

By the Bochner-Lichnerowicz formula, see e.g. [HUH], 

Dq = V*V + K, 

where TZ is the curvature endomorphism on V. Hence D 2 = V*V + 
B(x,D), where B(x,D) is a first order differential operator. Thus the 
hyperbolic Dirac-type equation in (|3*3"|) implies the wave equation 

(40) + V*V + B(x,D))u f (x,t) = 0. 

In the rest of subsection 3.1 we consider the local boundary condition 
associated to operator Pr 5 cf. (|28|) . Let u = u^(x,t), t G R + be the 
solution of 

(41) (id t + D)u(x,t) = inMxR + , 

Pr(u\dMxR + ) = /, u\t=o = 0, 

where / G C°°(dM x R + ,P r W). 

Since D r is a self-adjoint operator with V(D r ) C if 1 (M, V), the 
map U : f i— > it-* has an extension 

U:Hl{dM x [0, T], P r ^)^c7([0, T], # X (M, V)) n ^([0, T], L 2 (M, V)), 
when T > 0. 

Next we discuss the finite velocity of the wave propagation for the 
initial boundary value problem (jUJ. To describe it, let £ C dM be 
open. We define the domain of influence of E at time T as 

M(E, T) = {x G M | dist (x, E) < T}, 

where dist (x, y) is the distance in (M, g). Using local trivializations 
of the Dirac bundle and applying classical results on the finite velocity 
of the wave propagation in hyperbolic systems due to e.g. Wilcox |58j, 
we see that if supp (/) C E x [0,T], then 

(42) supp (u f (T)) C M(E,T). 

Finally, using the local boundary condition, we prove Lemma l2~5l 
Proof: Assume that we are given Cq{D) and P. Then Ap is determined 
by its graph {(?«,«) : u G Cq(D)} and visa versa. Since the induced 
bundle structure on dM determines Pr, Cq(D) determines Ap. □ 

3.2. Unique continuation. Our next goal is to supplement the finite 
velocity result by a unique continuation result of Tataru's type, see 
[53l l5i] for the pioneering works. To this end, for any open E C dM, 
we define the double cone of influence, 

K^ T = {(x, t) G M x R : \t - T\ + dist (x, E) < T} C M x R. 
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Theorem 3.1. Let £ C DM be open, u G C((0,2T), L 2 (M,V)) be a 
weak solution of the hyperbolic equation (id t + D)u = in M x (0, 2T) . 
Assume that 

(43) m|ex(o,2T) = 0. 
T/ien u(x,t) = for (x,t) G K(E,T). 

Proof: Assume first that u G C°°(M x (0,2T), V). Condition (jUJ) and 
the fact that w solves (f33l) imply that 

(44) 9jv«|ex(o,2T) = 0, 

where <9tv is the normal derivative. By (|i0|) . w satisfies the wave 
equation which, in a local trivialisation $ : it~ l {U) — > U x C d where 
U C M is open, may be written componentwise as 

(45) (<9 2 + D 2 )u = (<9 t 2 - A g )/u + Bi(a:, = 0. 

Here A g is the scalar Laplace operator on (M,g), I is identity matrix, 
and B\ is a first-order operator. We intend to apply Tataru's unique 
continuation result [HH1, in order to infer from (|i3|) . (jH)) , and (|4*5|l that 
it is zero in K(S,T). The difficulty is that the proof in [HB] deals only 
with a scalar wave equation. In the case of the scalar wave equation 
with time-independent coefficients a variant of the proof is given in |29[ 
Section 2.5]. This proof can be extended, word by word, to the vector 
wave equations of the form (<9 2 — A g )Iu + B\u = 0, where the leading 
order term is a scalar operator times the identity matrix. Therefore, 
g3J), JH, and (S3) imply that u = in K(E,T). 

Let now u G C((0, 2T), L 2 (M, V)) be a weak solution of the hyper- 
bolic Dirac-type equation (|33|) with the boundary condition (|4*3|) un- 
derstood in a weak sense. To prove the claim, we start by smoothing 
u with respect to time. Let tp{s) G C£°(— 1, 1), f ipds = 1. Then, 

f Is 

u £ (x,t) = / u(x,t - s)ip £ (s)ds, ip £ (s) = -ip(-), 

JR EE 

is in C°°((e, 2T — e), L 2 (M, V)) and solves there the Dirac-type wave 
equation {id t + D)u e = with u e |Ex(E,2r-e) = 0. As 

D 2 u £ = -d 2 u £ G C°°((e, 2T-e), L 2 (M, V)), 

the standard techniques of the elliptic regularity theory, e.g. (TJj, im- 
plies that u £ G C°°((e, 2T - e), H l (M, V)). Iterating this analysis we 
see that u £ G C°°{M x (e, 2T - e), V). By the first part of the proof, 

u £ = for {{x, t) : dist (x, E) + \t-T\<T- e}. 

As u £ — > u, when e — > 0, this implies that u = in K(S, T). □ 
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4. Inner products and controllability. 
4.1. Inner products. In this section we consider solutions u^(t) = 

u f (x,t) of (El-JHl, / e PC°°(dM x R + , W). We compute the norms 
of u J.(t) = n+u^(t) and u f _(t) = Yl_u f (t) using the boundary data, 
that is, the induced bundle structure (pTTj) on dM and the operator Ap. 
This type of formulae are called the Blagovestchenskii formula, due to 
their first appearance in |10| for the one-dimensional inverse problems. 

Using ((HJ) and ((T2j) we see that FD + DF = and thus U + D = 
DU_ and H_D = DH + . Thus, we can write the hyperbolic Dirac-type 
equation (|33|) as 

<«> M£ o))(:t:i)-° 

where = H + u(x, t) and U-(x,t) = U_u(x, t). 

Theorem 4.1. Let f,he PC°°(dM x E+, W 7 ) be given and u f , u h be 
the solutions of J^) -f^[j. Tften £ne response operator Ap, and the 
induced bundle structure on dM determine the inner products 

(47) ((u f + (t), <(s)>) and ((?/(t), M fe (s))) /ors,t>0. 

in particular, these data determine ((u^(t),u h (s))) for s, t > 0. 

Proof: We consider only the first inner product in (j%7J) as the second 
one can be evaluated in the same manner. The Stokes' formula for the 
Dirac-type operator D takes the form, 

(48) «D^»-«&-Ety»= / (a{N)<f>^)dA g , 

JdM 

where 4>,ip & if 1 (M, V) and oL4 g is the Riemannian volume on (dM, g). 
Denote by I(s,t) = ((U + u f (t),U + u h (s))) . Then, using ®, JI3), and 
(HE), we get 

(<9 2 - <9 t 2 )/M) = ((n+u'Cf), - ((u + dy(t),u h (s))) 

= ~((U +U f(t),D(Du h )(s))) + «Dn + ^(*),Dw h ( S ))) + 
-((DU +U f(t),Du h (s))) + ((flVW^"^))) 

= / « 7 (jv)n + u/(*), i>n + u fc (a)> + ( 7 (iv) J DnVW, n +M h ( s )» ovl, 
= i/ f( 7 (iv)n + A P /(t),n_a s (Ap/i( s )))- 

(49) -( 7 (iV)n_^(Ap/)(t),n + Ap/i(s))^ dA g . 

Since 7(0, t) = 1,(0, t) = l(s,0) = l t (s,0) = 0, we see that I(s,t) 
satisfies the one-dimensional wave equation lj%9~j) in the quarter plane 
s, t > with the known right-hand side and homogeneous initial and 
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boundary data. Thus, we can find I(s,t) = ((H + u^(t), U + u h (s))) for 
any s,i □ 

Remark 4.2. From the proof of Theorem 14.11 we obtain an explicit 
formula 

(50) ((IL±uf(T), u h {T))) = C ± [f, h]:= f [ X c(t, s)J^ s (f, h) dsdt. 

Here, the quadratic form C±[- , • ] is defined using the indicator function 
Xc is of the triangle C = {(t, s) G M 2 : T — t > \T — s\} and the 
quadratic form J± s (f, h), 

(51) 4 s (f,h) = l -f f( 7 (jv)n±Ap/(*),n F a.Ap/i( a )>- 

Z JdM \ 



-( 7 (iv)n T ^Ap/(t),n ± Ap/i( s ))J dA g . 

4.2. Global controllability. Let £ C dM be an open set and 

X P (S,T) = {«'(■ ,T) G L 2 (M,y) : / G PC °°(£ x (0,T), W)}, 

where u/(t) is the solution of (|33|) - (f3~ij) with the boundary condition 
associated with the operator P. In the following theorem, 

rad (M, S) = maxdist (x, £), rad (M) = maxdist (x, 9M). 

Theorem 4.3. i. Ze£ P be a boundary condition satisfying fTfy) and 
(Gjj andT > 2rad{M). Then the set X p {dM,T) is dense in L 2 {M,V) . 

2. Let Pr be the local boundary condition given by $28\) . £ C <9M be 
an open non-empty set, andT > 2rad(M, £). Then the set X r (E,T) = 
X Pr (S,T) 25 dense in L 2 {M,V). 

Proof: We start by proving claim 2. Let rj G L 2 (M, V) be orthogonal 
to X r (S,T). Consider the dual initial boundary value problem 

(52) (idt + D)v v (x,t) = inMxM, 

PrV v \dMxR = 0, V v \ t =T = V- 

Integrating by parts we see that, for any / G C£°(£ x [0,T], PrW), 

(53) = -i((uf(T),v v (T))) 

T r 

/ ((id t u f (t),v v (t)) - (u f {t),id t v v {t))) dVgdt 
((Du^f), - £>«„(*))) dV^d* 

(j(N)u f (x, t),v v (x, t))dA g {x)dt 
((j(N)P T u f (x,t),(I - P T )vr,(x,t)) + 



+ (j(N)(I - P T )u f (x,t),P T v v (x,t))) dA fl dt 



J M 
T 



./ftM 
T 
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(l(N)f(x, t),(I- P r )v v (x, t))dA g dt. 

JdM 

Here, the map 

(54) B : L 2 (M, V) ^ V'(dM x R, W), B( V ) = v v \ dMxR , 

is continuous. This may be shown by an approximation of r] G L 2 (M, V) 
by C£°(M, V) sections and applying the uniform boundedness princi- 
ple for the distributions. Since equation ijHHI) is valid for any / G 
C °°(S x [0,T],P T W), we see that (i"-PrK|sx(o,T) = 0. Together with 
the boundary condition of this makes it possible to apply Theo- 
rem l3.ll Hence v v vanishes in K^,t/2- However, the set contains 
an open neighborhood of Mx {^}. As also Prv v \dMxu = 0, this implies 
that v v = in M x R. In particular, 77 = f^(T) = proving claim 2. 

Claim 1. may be proven by the same arguments as claim 2. by 
taking P instead of Pr and using property (fl8|) . □ 

As can be seen in the above proof, for any t Q and / G C^°(dM x 
(0,t ), ^W), and v„ solving with 77 G L 2 (M, V), 

(55) «tt / (t ),^(*o)» = i /° / (7(N)f(x,t),Vr,(x,t))dA g dt, 

JO JdM 

where the right hand side is understood as a distribution pairing. 

4.3. Local controllability results. When P = Pr is the local pro- 
jection on dM, we define 

Xl(E,T) = {ny(-,T) G L 2 (M,V + ) : / G C °°(£ x (0,T), P r Vy}, 
A^(£,T) = {ILV(- ,T) G L 2 (M, V/_) : / G C °°(S x (0,T), iW)}. 

For an arbitrary set S 1 C M, we also denote 

L 2 (S, V) = {v G L 2 (M, V) : supp (v) C 5}. 

Theorem 4.4. Let Pr be the local projection on dM given by p%j) . 
T/zen, /or any T > and open S C <9M, £ne sets X±(£,T) are dense 
in the subspaces Tl±L 2 (M(Y*, T), V) C L 2 (M,V), correspondingly. 

Proof: We will prove this result for X£(E, T), the claim for X^(E,T) 
may be obtained in the same manner. 
Let T) G n+L 2 (M(S,T),1/) satisfy 

((77, II + </(T)» = for any / G C °°(S x (0,T), P r W). 

Consider the initial boundary value problem 

(56) {id t + D)v v {x,t) = inMxK, 

PrVnldMxR — 0, 
n + n^|i =T = 77, n_n^|t =T = 0. 
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Then v v G C(R, L 2 (M, V)) and, as noted in the proof of Theorem Q| 
it has a well defined trace v v \qmxr G V(dM x R). Same considerations 
as in (f53j) show that 

0=/ / (7{N)f{x,t),{I-P r )v v {x,t))dA g dt, 

JO JdM 

for any / G C °°(Ex [0, T], iW). Therefore, (J-P r )^|sx(o,T) = 0. To- 
gether with the boundary condition in (f56|) . this imply that ^|s x (o,T) — 
0, i.e. H±v ti \y:x(o ,t) — 0. Using (fi6|) and the fact that U.^v v \ t= T — 0, we 
see that n_v^(T + s) = — Il_t> r? (T — s) and Il + ^(T + s) = U + v v (T — s). 
Therefore, 

(57) n ± ^| Sx(Ti 2T) = and supp (^|ex(o,2T)) C S x {T}. 

We intend to show that v v \qmx{o,2T) = 0. Due to the boundary condi- 
tion in (|5H this is equivalent to 

(58) / [ (j(N)f(x,t), v v (x,t)) dA g dt = 



for any / G C °°(S x (0, 2T), P T W). In view of (JSZJ, it is enough to 
prove that 

(59) / / (-y(N)f e (x,t), v v (x,t)) dA g dt 

JT-2e Jy, 

= i((u fs (T + 2e), ^(T + 2e))) — ► 0, when e ->• 0. 

Here f E (-,t) = Xe(t — T)f(-,t) and Xe( s ) is a smooth cut-off function, 
X {s) = 1 for s G (-1, 1), supp(x) C (-2, 2) with ^(s) = x(s/e). Note 
that the equality in (f59|) follows from (|55|) . 

Assume, in the beginning, that f E (x,t) = Xe(t — T)Fg(x), where 
Fa = F\q M ,F G C°°(M, V). We can represent the solution as 

u f c{x,t)= X e(t-T)F{x)+uj £ {x,t), 

where u £ is the solution to 

(id t + D)u £ = -{id t + D) (XsF) = --x 1 ((t - T)/e) F - Xe (t - T) DF, 

e 

Pr^e\dMxR = 0, U £ \t=T-2e = 0. 

We represent F and DF as L 2 (M)-converging series of eigenfunctions 
4> k of D T , 



oo 

L fe<pfe- 



F = J2^k, DF = J2^ 

fc=i k=i 

Then 

oo 

(60) u £ {T + t) = Y, H(t) + w 2 k (t)) <p k = u]{T + t) + u 2 £ (T + t), 

k=l 
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wl(t) = - f e lX ^ X \s/e)ds, w 2 k (t) = -b k f e iX ^ X (s/s)ds. 

IS Jt-2e JT-2e 

It then follows from (jfiUl) that 

(61) \\ul(T + 2e)\\ L 2 < d||F|| L2 , \\u 2 £ (T + 2e)\\ Ifi < C 2 e\\DF\\ L 2. 
What is more, introducing 

ul' N (T + t)= wlWfa, 

\\ k \>N 

we see that 

(62) lim || £ j e 1 ' JV (T + 2e)|| L2 = 0, 

uniformly with respect to e G (0, Sq). To estimate u)\ — uj\' n , use 
integration by parts in the equation for w\{t) in (|60f to get 

(63) K(2e)| < C 3 N\a k \e when |A fc | < N. 

Combining f^-JEH) we see that \\uo l £ (T + 2e)\\ L i -> when e -> 0. 
This and the estimate for (T + 2e) in (RUT) imply the desired result 
©for/ 6 (x,t)=x e (t-r)F 8 (a;). 

Returning to a general / G C£°(£ x (0, 2T), PpW), it remains to 
estimate uj £ {T + 2e) for / e (x, t) — (t — T)x £ (t — T)f(x, t) with smooth 
f(x,t). Arguments similar to those leading to the estimate of oj 2 {T + 
2s) in (JHH) show that ||c<; e (T + 2e)||x,a — > when e — > 0. This proves 
and, therefore, Theorem 14.41 □ 

We complete this section by describing a procedure to find rad (M, E). 

Corollary 4.5. For any open S C <9M, the Cauchy data set C (D) 
and the induced bundle structure on dM determine rad(M, E). 

Proof: By the local controllability result, Theorem 14.41 we see that, if 
Ti > rad (M, E), then the following claim (A) is valid: 

(A) For any T 2 > and / G C^{dM x R+, P r VF) there are hj G 
C °°(S x (0, Ti), iW), j = 1, 2, . . . such that 

(64) lim ||n + (« ft i (TO -« / (T 2 ))|| z2 = 0, 

lim ||n_(tt fc 7(T0-« / (T 2 ))|| i3 = 0. 

On the other hand, if T < rad (M, E), we see that for T 2 > rad (M, E) 
there is / such that supp (u^(T 2 )) (jL M(E,T). Therefore, claim (A) is 
not valid. 

As by Theorem 14.11 the norms on the left-hand sides of (|6*4*jl can be 
evaluated using the boundary data, we find rad (M, E) by taking the 
infimum of all Tj > such that (A) is valid. □ 
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5. Generalized boundary sources 

5.1. The wave norm. We define the wave operators 

U T : PC£°(dM x (0, T), W) — > L 2 (M, V), U T f = u f (T), 

where T = 3 rad (M) and u f (t) = u f (x, t) is the solution of (1331) -CT. 
By means of the wave operators, we define a semi-norm, or the wave 
semi-norm on the space of boundary sources as 

(65) \\f\W = \\U T f\\ LHMy) . 

By Theorem 14.11 the knowledge of the operator Ap and the induced 
bundle structure (pTT| on dM enables us to calculate explicitly the semi- 
norm JHSl) of any / G PC^(dM x (0, T), W). 

To complete the space of boundary sources, we define the equivalence 
relation ~ on this space by setting 

/ ~ h if and only if u f (T) = u h {T). 

Further, we define the space B as 

B = PC™{dM x (0, T), W)/ ~ . 

Then (f65|) becomes a norm on B. Finally, we complete B with respect 
to the norm (JfiHj) . This space, denoted by T = B, consists of the 
sequences of sources which are the Cauchy sequences with respect to 
the norm (|65|) . We denote these sequences by /, 

/ fjePC?{dMx(Q,T),W)- 

Also, we denote J : PC£°(dM x (0, T), W) — > T the factor-projection 
that maps / to its equivalence class in B. 

The sources / G T are called generalized sources in the sequel. The 
corresponding waves are denoted by 

(66) Wf(t) = lim u fj (t) fort>T, 

where the convergence in the right-hand side is the L 2 (M, V)— conv- 
ergence. By the global controllability, Theorem 14.31 and the definition 
(jfiHJ) of the norm, the map / h- > Wf(t) is an isometry, T — > L 2 (M, V). 
We extend the notation Wf(t) for all t G M by defining 

(67) Wf(t) = exp(i(t - T)D P )Wf(T) for t G E. 

Due to this formula we call Wf(t) a generalized wave. 

Remark 5.1. The choice of the parameter T = 3 rad (M) in the above 
construction can be replaced with any fixed T > 2 rad (M). 

We note that the above construction of the space of generalized 
sources is well-known in PDE-control, e.g. [4(H 150]. 



INVERSE PROBLEMS FOR DIRAC OPERATORS 21 

For t G R and any boundary source / G J 7 , we can find if the 
condition 

(68) 3W/(t) G L 2 (M, V) 

is satisfied for a given I G Z + . Indeed, for / = 1, condition (|68|) for 
/ — (fj)fLi G is equivalent to the existence of h E J 7 such that 

(e9) h ,»zMM)_ m(f)li , (1 , fl 




= lim||W| — ^-A ](/)|U,,„., , II. 

t > T. Here 7^ : f(x,t) i— > /(x, t — 5) is the time-shift naturally 
extended from PC^°{dM x (0,T),W) to J". Note that the validity 
of (|69| is independent of value of t > T as the norms ||W/(t)||ia are 
independent of t. Using the induced bundle structure on DM and 
the response operator Ap we can verify if a given generalized source 
/ — (fj)'jLi an d a trial h = {hj)JL 1 satisfy condition (|69j) . In this case 
we denote 

h = dj and feV(d t ). 

Having defined d k ~ l f we define d k f, if it exists, by induction. The 
spaces of generalized sources which have k G Z + time derivatives will 
be denoted by T k = V(d$) with J 700 = f\ eZ+ T k . Clearly, 

dtJf = Aft), f G PC™(dM x (0, T), W), 

that is, d t is an extension of the time derivative f ^ ft defined in the 
usual way for / G PC^{dM x (0,T), W). Note that by (EZJ) and @J$, 

(70) / G V(d k ), if and only if Wf(t) G V(D k P ), t G R. 
In addition, 

(71) W{d k J){t) = i k D k P Wf{t) for / G 

5.2. Boundary sources corresponding to eigenfunctions. When 
/ G T x we can evaluate, using Theorem 14. 1| the inner products 

((D P Wf(T), W/(T)» = -i((W(d t f)(T), W/(T)», 
({D P Wf(T), D P Wf(T))) = ((W(dJ)(T), W(dJ)(T))), 

as well as the inner products ((Wf(T), Wf(T))). On the other hand, 
by the Hilbert-Courant min-max principle the eigenvalues [Aj] 2 of Dp 
can be found as 

(72) [Xj] 2 = sup inf (( D ^ D ^\ 

where Uj are in V(Dp). This makes it possible to prove the following 
result: 
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Proposition 5.2. The induced bundle structure JUJj on dM and the 

response operator A p determine the eigenvalues Xj of Dp and the gen- 
eralized sources hj G J 700 , j = 1,2,..., such that <fij = Whj(T) 
are orthonormal eigenf unctions of Dp that form a complete basis of 
L 2 {M,V). 

Proof: Using (|72^l . we can find all the eigenvalues of Dp. Indeed, by 
the global controllability, Theorem 14.31 and relation (|7T|) . W/(T) runs 
over V(Dp) when / runs over T x . At the same time, the minimizers in 
(j?2*J) are those generalized sources fj G T 1 that the waves Uj = Wfj (T) 
are the corresponding orthonormal eigenvectors of Dp. Furthermore, 
since we can evaluate ({D P Uj, Uj)), we can decompose the eigenspace 
of Dp, which corresponds to an eigenvalue [X,] 2 , into a direct sum of 
the eigenspaces of Dp which corresponds to the eigenvalues Xj and 
—Xj. Hence, we can find hj G T x such that 4>j — Whj(T) are the 
orthonormal eigenfunctions of Dp. Clearly, hj G T°° . □ 

Applying Theorem 14.11 to W/(i) and Whj(T), we can evaluate the 
Fourier coefficients Wj(t) = ((Wf{t), fa)), j = 1, 2, . . . , of W/(t). Us- 
ing Wj(t), we can compute the norms ||(1 + \D P \) s Wf{t) \\l 2 (m), s >0. 
We define the spaces T s as the spaces of those generalized sources 
feF for which Wf(t) G V(\D P \ 8 ) and define 

oo 

ll/V = || (i + \D P \ywf(t)\\ LHM) = + \Xj\r\wj(t)\ 2 , 

3=1 

where the right-hand side is independent of t G R. 

6. Index formulae using boundary data 

In this section we use boundary data to verify if the decomposition 
condition (|24j) is valid and give formulae for the Fredholm index of the 
operator Dp appearing in decomposition (|2*51) . 

6.1. Decomposition of D P . Let us consider first the Dirac operator 
D r associated with Pp. Using properties (JBJ) of the chirality operator, 
we can see that, if u G V(Dp) and U + u G V(Dr), then u\qm = 0. 
Therefore, the decomposition condition (|2*ljl is not satisfied for the 
Dirac-type operator Dp. However, there is a large number of important 
examples of projectors P when the decomposition condition ((2H) is 
valid for the corresponding Dirac-type operator Dp. Some of these are 
discussed below. 

6.2. Decomposition of T. Assume that we are given the operator Ap 
with a possibly non-local boundary operator P and the induced bundle 
structure on dM. Then the formula (jHOjl determines the non-negative 
sesquilinear forms £+[■,■] and £_[•,•] on PC^°(dM x (0, T),W) such 
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that 

C ± (f,h) = ((U ± uf(T),u h (T))). 

Observe that the quadratic form C = C + + £_ defines a semi-norm 
y/C[f, J] that coincides, on PC™(dM x (0, T),W), with the semi-norm 
defined by (|65|) 

(73) C[f,f} = C + [f,f}+C4f,f] = \\f\\% 

Therefore, C + and £_ can be extended to T . 

In turn, these forms give rise to two bounded non-negative operators 
B± on T so that 

The fact that U + L 2 (M,V) © Tl^L 2 (M,V) = L 2 (M,V) implies that 
B± are orthoprojections in T satisfying B + © £>_ = I. Introducing an 
operator S, 

Sf = Wf{0), feF, 

we have 

(74) SB J = U ± Sf. 



6.3. Index formulae. At first we prove Theorem 12.61 
Proof: Using the operator S, equation (|7lT) may be written as 

(75) Dp = S(-td t )S-\ 

Together with (|74| . this equation imply that the condition Yi±D(Dp) C 
V(Dp) is satisfied if and only if B±D(d t ) C V(d t ). This proves claim 
1. 

If decomposition (|2^| is valid, so that Dp can be decomposed as in 
(|25|) . formula (fT5|) implies that d t can be decomposed as 

dt= {d+ \ ) -B+V^^B-V^^B-FQB+F, 

where df are the restrictions of d t on B±U(d t ), correspondingly. As 
Dp are Fredholm operators, we conclude from (ffUl and (fT5|) that df 
are Fredholm operators and Ind (df) = Ind (Dp). This proves claim 
2. ' □ 

Corollary 6.1. Assume that condition p^p is valid for Dp. Then 

Ind (Dp~) = sgn (\B + - S_]| Ke r(s t )) • 

Proof: As FD + DF = 0, we have F : Kev(D P ) Ker(D P ) so that 
U± : Ker(Dp) -> Ker(D P ). Thus, Ind (<9 t + ) is equal to 

dimKer (9t| B+ ^i) - dim Ker (d t \ B _ T i) = sgn ([£+ - 5_]| Ker(at) ) , 

where we use the fact that _B± are orthoprojectors in T . □ 



24 



YAROSLAV KURYLEV AND MATTI LASSAS 



By Proposition \5.2\ using the boundary data it is possible to find 
generalized sources hj = (h^)^ =1 G T°° such that Whj(T) = <f>j, with 
{4>j}'jL l being an orthonormal basis of the eigenfunctions of Dp. Let 
4>j, j = 1, . . . , v, where v = dimKer (-Dp), be the eigenfunctions span- 
ning Ker (Dp). Then, for any to > T, 

lim Ap(^)|aA/x(t ,oo) = l(t) <fij(x)\dM, j = l,...,v, 
with l(t) = 1. Thus, Ap determines the subspace 

Af = span^U : j = 1, ■ ■ ■ ,u} C (1 - P)C°°(dM, W). 

As F : Af — > Af and F = Il + — Corollary 16. II implies the following 
result: 

Lemma 6.2. Assume that Dp satisfies condition $%\) . Then 

Ind(D+) = sgn{F\ M ) , 

where the right-hand side can be evaluated in terms of the induced bun- 
dle structure on dM and the response operator Ap. 

Example 1, continued. It follows from the definitions of T and the 
relative and absolute boundary conditions, that the Dirac operators 
(d + 5) r and (d + 5) a can be decomposed into (d + 5)f and (d + 5)^. 
Clearly, Ker[(c? + 5)+] consists of the relative harmonic forms of the 
even order and Ker[(<i + 5)~] of relative harmonic forms of the odd 
order. Similarly, Ker[(d + <5)+] and Ker[(<i + S)~] consists of absolute 
harmonic forms of an even and odd order, correspondingly. Therefore, 
Ind [(d + <5)+] is the Euler characteristic x(M) of M. By the Poincare 
duality, Ind [(d + 5)+] = ±Ind [(d + 5)£] depending on n = dim(M). 
Thus, Theorem 12.61 gives a representation of the Euler characteristic 
x(M) in terms of the Cauchy data set Co(D) of the hyperbolic form- 
Dirac equation. 

Example 2, continued. Similarly to example 1, the index of the ex- 
tended Maxwell operator i(d — S a ) with the electric boundary condition 
is the Euler characteristic x{M). 

Example 3, continued. The generalized signature operator on even 
dimensional manifold with boundary satisfies decomposition condition 
(f2i|) . This was first considered in the famous paper [l] (for further 
developments in this area see e.g. [121 CHI SU EH])- Theorem 12.61 
provides a possibility to find the index of the signature operator from 
the Cauchy data of its Dirac operator. 

7. Reconstruction of the manifold 

In this section we will show how to reconstruct the manifold, M and 
metric, g from the response operator Ap. To this end, we will show that 
the boundary data determine the set of the boundary distance functions 
that determine (M, g) up to an isometry. 
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7.1 Controlling supports of generalized waves. We start by fixing 
certain notations. Let 

(76) Ti = T + 4 sup {rad (M, E) : E ^ 0}. 

Then Ti is determined by the boundary data and Ti — T > diam (M) . 

Let Ej C <9M be non-empty open disjoint sets, 1 < j < J and r~ 
and be positive times with 

< tJ < r/ < Ti - T, 1 < j < J. 

Let 5 = SdEj, Tj - , r J + }^ =1 ) C M be the intersection of slices, 

j 

(77) 5=fl(M(E,,r;)\M(E„rr)). 

i=i 

Our first goal is to find, using the boundary data, whether the set 
S contains an open ball or not. To this end, we give the following 
definition. 

Definition 7.1. The set Z = Z({Y,j,rj~ ,Tj~}j =1 ) is the set of the gen- 
eralized sources / G that generate the waves Wf (t) with 

(1) supp (W/(Ti)) C M(Ejf, r+) for all j, 1 < j < J, 

(2) W/(Ti) = in M(Ej, r") for all j, 1 < j < J. 

Theorem 7.2. Given the induced bundle structure on dM and the re- 
sponse operator A?, we can determine whether a given boundary source 
f G T°° is in Z or not. 

Proof: Let / = (fk)kLo G be a generalized source. Consider 
first the question whether supp (W/(Ti)) C M(Ej, Tj~) or equivalently, 

supp(n + W/(Ti)) C M(E J ,r+) and supp (LW/^)) C M(E i ,r/). 
By Theorem l4.4l this is equivalent to the existence of generalized sources, 

h + = (hj)Zo and h- = (hi)Z , hf e C^(E j x (0, r+), P r W), 
such that 

(78) lim \\Yi + {u^{T{)-u^{r+))\\ L , {My) = ^ 

(79) Jrn^ Hn-^CTO - u^(r+))|U W ) = 0. 

By Theorem 14. 1| we can evaluate the left-hand sides of (f78|) and (f79| 
using the induced bundle structure (j3TJ) on dM and the operator Ar. 
Thus, for a given /, we can verify if condition 1. is satisfied. 
Next, consider condition 2. As 

(id t + D)Wf = in M x R, P r (w/| 9M ) = 0, 

condition S. implies, due to the finite propagation speed, that 

Wf = in = {(x, t) G M x R : dist (x, E,-) + |t - Ti| < r~}, 
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for all j = 1, ... , J. In particular, 

(80) W/i E . x(Tl _ T - jTl+rJ) = lim Ar/,-| s . x ( Tl _ T -, Tl+ rr) = 0. 

On the other hand, if (jKUj) takes place, then Theorem 13.11 implies that 
Wf = in Kj. As condition (j80|) can be verified given the induced 
bundle structure on dM and the response operator Ar, we can verify 
condition 2. □ 

Theorem 7.3. Let S and Z be as above. The following alternative 
holds: 

1. If S contains an open ball, then dim(Z) = oo, 

2. If S does not contain an open ball, then Z = {0}. 

Proof: 1. Let B C S be an open ball and f ip G C™(B, V). By the 
global controllability, Theorem 14 .3\ there is / G J 700 such that 

W/(Tx) = p. 

As Q°CB,V) C V(D^), then / G Z. As G C^(B,V) \ {0} is 
arbitrary, this implies that dim(Z) = oo. 

2. Assume S does not contain an open ball and suppose that there 
is 7^ / G Z. Then, by conditions 1., 2. in Definition l7.lt 

J 

supp(W/(T 1 )) C f| M(S„r+) = 5 + , 

Wf{x,T l ) = 0, for x G U/ =1 M(S J; rr) = 5-. 

Thus supp(W/(Ti)) C S + \S~ . As 5 does not contain a ball, then S + \ 
S~ is nowhere dense. Since W/(Ti) is smooth, it vanishes everywhere 
in M. This is a contradiction with the assumption that / 7^ in T . □ 

We are now ready to construct the set of the boundary distance func- 
tions. For each x G M, the corresponding boundary distance function, 
r x G C(dM) is given by 

r x : dM -»• R + , r x {z) = dist (x, z), z G 9M. 

The boundary distance functions define the boundary distance map 1Z : 
M — > C(dM), TZ(x) = r x , which is continuous and injective (see [HH 
ES|). Denote by 

K{M) = {r x G C(dM) : x G M} 

the image of IZ. It is known (see jHHEll) that, given the set TZ(M) C 
C(dM), we can endow it, in a constructive way, with a differentiable 
structure and a metric tensor g, so that (TZ(M),g) becomes a manifold 
that is isometric to (M,g), 

(TZ(M),g) = (M,g). 
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We complete the reconstruction of (M,g) with the following result: 

Theorem 7.4. Let the induced bundle structure on dM and the 
response operator Ap be given. Then, for any h G C(dM), we can find 
whether h G 1Z(M) or not. 

Proof: Consider a function h G C(dM). We have that h G K{M) if 
and only if there is x G M such that h = r x . Observe that if such x 
exists then, for any e > 0, we have h(x) — e < d(y, z) < h(x) + e for 
any z G dM and y G B(x,e), where B(x,e) C M is x-centered ball of 
radius e. Therefore, for any J G Z+, any set of points Zj G dM, and 
non-negative numbers rj, T+ satisfying rj < h(zj) < Tj~, j — 1, . . . , J, 
the set 

(81) S = S({Hj, r~, rj~}j =1 ) C M contains an open ball, 

if Ej C <9M are sufficiently small neighborhoods of of Zj. 

On the other hand, assume that for all J, Zj G dM, r~, satisfying 
rj" < h(zj) < Tj~, and sufficiently small E^, condition (|8T|) is satisfied. 
Then there is a point x G M such that 

\d(x, Zj) - h(zj)\ < r+ - rr + diam (Ej), j = 1, 2, . . . , J. 

Thus, letting r~ — > — > and, diam (Ej) — > 0, we see that 
there is xj G M such that = h(zj) for all j = 1, 2, . . . , J. 

Denote by {zj}JL 1 a countable dense set on dM, and let xj, J G Z + 
be the above defined points corresponding to the first J points Zj, j = 
1, . . . , J, of this set. This gives us a sequence (xj)f =1 . As M is compact, 
this sequence has a limit point x G M. By density of {zj}°^ 1 , h = r x . 

By Theorem I7.3| condition (|8*TT) is equivalent to the condition 

dim(Z({E„T J -,r+}/ =1 ))^0. 
By Theorem \7.2\ this can be verified using the boundary data. □ 
Combining the above theorem with Lemma 12.51 we obtain the fol- 
lowing result: 

Corollary 7.5. The induced bundle structure on dM and the 
Cauchy data set Cq(D) determine the Riemannian manifold (M,g) up 
to an isometry. 

8. Focusing sequences. Proof of Theorem 12.91 

8.1 Waves and delta-distributions. In the previous section it was 
shown that, using the induced bundle structure on dM and the re- 
sponse operator Ar, we can control the supports of the waves W/(t). In 
this section the goal is to construct a sequence of sources, (fi)^ C 
such that, when I — > oo, the corresponding waves Wfi(t) concentrate, 
at t = Ti, where T\ is defined by (fTfil) . at a single point y G M mt . 
These waves turn to be useful for reconstructing the bundle V and the 
Dirac-type operator D. 
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In the following, let S y denote the Dirac delta-distribution at y G 
M int , i.e., 

/ 6 y (x)<f>{x)dV g {x) = cf>(y), for G C °°(M). 

Since the Riemannian manifold (M, 51) is already found, for any y G 
M mt we can choose a sequence C dM, < < r+ < 7\ — T, so 

that, for 5i = 5({S j7 , r+, j^f J), 

(82) S l+1 cS h = 

1=1 

Let = Z({Eji,rJ l ,T^}j^]) be the corresponding set of the general- 
ized sources defined in Definition 17. 11 

Definition 8.1. Let S h I = 1,2,..., satisfy (JH2J) . We call the se- 
quence (fi)™i, with fi G Zi, a focusing sequence of generalized sources 
of the order s G R+ (for brevity, a focusing sequence), if there is a 
distribution- valued section 7^ A y G V(M, V) supported at the point 
y G M in * such that 

lim(<W/,(7\)^>> = ((A y ,r])), for all V G Pd^rl 8 )- 

Lemma 8.2. Assume that the induced bundle structure on dM 
and the response operator Ar be given. Then, for any s G M + and 
any sequence of generalized sources, (fi)fZx c •^ r °°j w,e can determine if 
(fi)™i « s a focusing sequence of the order s or not. 

Proof: By Theorem 14.31 any r\ G V{\Dy\ s ) can be represented as r\ = 
m(Tx), /i G T s . Thus, if (/,)£ x is a focusing sequence of the order s, 
then the limit, 

(83) ((Ay, Wh(Ti))) = ]jm((wMTi), Wh(Ti))), 

I— >oo 

exists for all h G T s . By Theorem 14. 1| the existence of this limit can 
be verified if we are given the induced bundle structure on dM and the 
operator Ap. 

Conversely, assume that we have a sequence (/i)£i> with f\ G Z\ such 
that the limit (|83|) exists for all ) k E T s and that this limit is non-zero 
for some h G T s . Then, by the principle of the uniform boundedness, 
the mappings 

ri^((Wfi(Ti),v)), leZ+, 
form a uniformly bounded family in the dual of V(\D r \ s ). By the 
Banach-Alaoglu theorem, we find a weak* -convergent subsequence in 
this space with the limit 

A y := KmW/iCZi) G (v(\D T \ 8 )] . 
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This gives us the required non-zero distribution- valued section in V(M, V). 
As fi£ Z h then supp(A y ) = {y}. □ 
Next we consider the representation of A y in a local trivialisation 
near y. Since supp(A y ) = {y}, A y is a finite linear combination, with 
coefficients from 7r _1 (?/), of the Dirac delta-distribution at y and its 
derivatives. The role of the smoothness index s is just to select the 
order of this distribution, as is seen in the following result. 

Lemma 8.3. Let A y , y G M mt be a limit of a focusing sequence of the 
order s, | < s < . Then A y is of the form 

A y (x) = X5 y , 

where A G n~ l (y). Furthermore, for any X G n^ 1 {y) there is a focusing 
sequence suc ^ that the waves Wfi(Ti) — > X5 y in V(M, V) when 

I — > oo. 

Proof: As for s > 

# c s omp (M, V) C V(\D r \ s ) C H S {M, V), 

the embedding theorems imply that, for | < s < the delta- 

distributions XS y , X G 7r _1 (y), are in (T>(\D r \ s ))' while the derivatives 
of the delta-distribution are not in (V(\D r \ s ))' . 

To prove the assertion we have only to verify that a desired focusing 
sequence exists for any A G n^^y). This can be done by an approxi- 
mation of XS y by smooth sections with supports converging to y and 
an application of the global controllability Theorem 14.31 □ 
8.2. Construction of the bundle structure. 

Lemma 8.4. Assume that we are given the induced bundle structure 
pI7\) on dM and the response operator Ap. Then, for any yo G M mt , 
we can find integers k(l), . . . , k(d) and a neighborhood U C M mt of yo 
such that the eigenf unctions 4>k(i)(y), ■ ■ • > 4>k(d)(y) of form a basis in 
^^iy), V e U. Moreover, for any h G J 700 , we can find the coefficients 
cti(y,t) of the decomposition 

d 

(84) Wh(y, t) = J2 a ^ *)&(i)(3/)> yeU,teR. 

i=i 

Proof: Let Ui C M mt be a sufficiently small neighborhood of a point 
y G M mt so that there is a trivialisation $ : 7T —1 (Z7i) — > U-y x C d . 
For all y G U we choose d families of focusing sequences, (fl' y )kLi, 
j = 1, . . . , d, such that Wfl' y (T) converge to some distributions X J y 8 y , 
A^e - ■(!/). ytl\. 

Using the operator A r and the induced bundle structure on dM we 
can evaluate, by Theorem 14. 1[ the inner products 

(85) lim((Wh(T),Wfl> v (T))) = ((Wh(T), = (Wh(y,T),\ 3 y ) y , 
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where j = 1, . . . , d. The set {Wh{T) : h G J 700 } contains C °°(M, V). 
Therefore, {Wh(y, T) : h G = 7r _1 (y) and it is possible to verify, 
using JHSJ, whether the set {A^, j = 1, . . . , d} is a basis in 7r _1 (y) for 
y lying in some neighborhood U C t/i of ?/o- Moreover, we can verify 
whether A^, considered as functions of y, define C^-smooth sections in 
U. 

On the other hand, Lemma IP1 guarantees the existence of focusing 
sequences such that the corresponding X J y are linearly indepen- 

dent and smooth in U. 

For any h G J 700 , i El, consider the map K t :T°°^ C°°{U, C d ), 

K t h(y) = ((Wh{y, t),X l y ) y ), . . . , (Wh(y, t), \ d y ) y ), y e U. 

By Proposition l5.2[ we can find the generalized sources hk, k = 1, 2, . . . 
such that Whk(T) = Thus, for any o? numbers fc(l), . . . , k(d) G Z+, 
we can evaluate the matrix 

:= [ e jv(2/)]i=D e iAv) = (4>Hi)(y), Mj) y , y e Z7i- 

As for any y £ A^ mt the vectors {4>k(yo)}kLi span 7r _1 (|/o), it is possible 
to select k(l), . . . , k(d) so that the system $ of d eigenfunctions, $ = 
• • • ,4>k{d)) is linearly independent at yo- As £(?/) is smooth, we 
can choose U C U± to be a neighborhood of yo where is invertible. 

Let h G Then the coefficients cti(y,t) in can be found by 
computing the vector (ai(y, t), . . . , a d (y, t)) = £(y)~ 1 K t h(y). □ 

Next we consider an eigenfunction basis near the boundary. Due 
to the boundary condition Pr(4>k\dAi) — 0, the eigenfunctions <fik do 
not span 7r _1 (|/o), yo £ dM. Because of this, we will use near dM a 
combination of the eigenfunctions 0^ of D r>+ and (f)J of D T -, where 
D r + = D r and -Dr,- = Dp with P = (I — Pp) (compare with [2] where 
a coordinate system near dM was constructed using a combination of 
the Dirichlet and Neumann eigenfunctions). We denote the spaces of 
the generalized sources corresponding to the operators Dr,± by T± and 
F± and by W± be propagators (ffifil) for Dp±. We note that our previ- 
ous results concerning the unique continuation, controllability, and the 
construction of focusing sequences can be proven using the boundary 
condition I — P T — |( J — r) instead of Pr = \{I + T). 

To simultaneously consider the operators D r + an d -Dr,-, we need 
the following modification of Theorem 14.11 

Lemma 8.5. Assume that we are given the induced bundle structure on 
dM and the response operator A r . Then for f + G C£°(dM xR + , P r W) 
and f- G C oo (9Mxt + , (l—Pp)W) we can evaluate the inner products 

((u f +{t),u f -(s)}}, fort,s>0. 

Here is the solution of with P = Pr and is the 

solution of CT -OT with P = I - P r . 
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Proof: Let I(s,t) = ((u^+(t), i^-(s))). Similarly to the proof of Theo- 
rem HH1 we use integration by parts to see that 

(86) (d s + d t )I(s,t) =i((Du f+ (t),u f -(s)))-i((u f+ (t),Du f -(s))) 

= i[ ( 7 (JV)A+/ + (i) I Af/_( s ))^. 

JdM 

Here Ap are the response operators for D r ±. By Lemma l2~5l the op- 
erator Ar = Ap" and the induced bundle structure determine Ap. As 
1(0, t) = I(s, 0) = 0, we can determine I(s, t) by solving the differential 
equation (|86|) along the characteristics t — s — c. □ 

Let us denote by tpk °f the eigenfunctions of D r + and -Dr,-, renum- 
bered so that ip2k-i = 0^ and ip2k = 0& > k > 1. To generalize Lemma 
18.41 we need the following result: 

Lemma 8.6. V(Df? + ) + £> (£>{?_) = C°°(M, V). 

Proof: For 2 G dM let be a sufficiently small ball B = B(z,r) C M 
such that the boundary normal coordinates and a local trivialization 
$b : 7r _1 (i?) — > B x C d are well defined in £?. It suffices to show that 
any u G C°°(M, V) with supp(w) C B(z,r/2) can be written as 

(87) « = + where u + G £>(£>;?+), w_ G V(D™_). 

To prove JH3), let u J G C°°(<9M, W 7 ), j = 0, 1, 2, . . . be supported on 
B(z,r/2) fl dM. By using Seeley's extension theorem [HI] in B, there 
is w G C°°(M, V) supported on B(z,r/2) such that 

difiilsM = u 3 . 

Moreover, using representations (J2J), (fl6j) of the Dirac operator in the 
local trivialization we see that 

(88) D%U = 7(#V + A 3 [u\u\ . . .y- 1 ), 

where A, is a differential operator. Let G C^(dM, W), j = 0, 1, 2, . . . , 
be supported in B(z, r/2) HdM. Solving recurrently equations (|8*%|) we 
see that there is u G C°°(M, V) such that D^uIqm = Wj. 
Let v G C°°(M, V"). Then v G £>(£>^ ± ) if and only if 

(89) Pr,±(i}V)U = 0, j = 0,1,2,.... 

Thus, if it G C°°(M, V) supported in B(z, r/2) is of the form JEH, we 
have by ((HJ) that 

£>VU = (1 - PrX^M =: 

D J u-\aM = Pr{D 3 u\dM) =■ wj. 

On other hand, let u G C°°(M, V) and G C°°(M, V) be func- 

tions such that D^u±\qm = wf. Then u± G £>(-D£° ± ) and (j%7jl is 
satisfied. □ 
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Lemma 8.7. Assume that we are given the induced bundle structure 
( I Hl\) on dM and the response operator Ap. Then, for any yo G M , we 
can find integers k(l), . . . , k(d) and a neighborhood U C M of yo such 
that the eigenf unctions ipk{i){y)t ■ ■ ■ > V'fe(d) (2/) °f Dr,+ or Dr,~ form a 
basis of n~ l (y), y G U . Moreover, for any h G J-± , we can find 
coefficients af{y,t) such that 

d 

(90) W±h ± (y,t) = ^actfaWwiy), yeU,teR. 

i=i 

Proof: First we show that, for any y G dM, the space 7T _1 (yo) is 
spanned by {^(j/o)}^=i U {fa (yo)}k=v To this end, let u G C°°(M, V). 
By Lemma ESI there are u + G V(D^ + ) and w_ G V(D^_) such that 
w = u + + tt_. Consider the eigenfunction expansions 

oo 

(91) M± = ^a J ± 0± 

fc=i 

that converge in X>(-Df±). By elliptic regularity in Theorem 12.11 the 
series (|9"T|) converges in C°°(M, V). In particular, any jj, G 7r _1 (y ) can 
be written in the form // = ajfa(y ) + EfcLi ajfa(yo)- 

Now, let |/o G 9M. Arguing as in the proof of Lemma [H31 but using 
Lemma instead of Theorem I4.1| we find k(l), . . . , fc(rf) G Z + and 
neighborhoods C/ of yo such that ipk(i)(y), ■ ■ ■ , i>k(d)(y) form a smooth 
basis in 7r -1 ({7 fl M int ). Moreover, we can find, using the boundary 
data, the coefficients af(y,t) in <$ty for any y G U H M int . Then, 
by Lemma IKTHl ipk(i)(y), ■ ■ ■ ,i>k(d)(y) form a smooth basis in 7r — 1 (C7) if 
and only if the coefficients af(y,t) may be continued smoothly on the 
whole U for any ft, G JF^ 3 . This can be verified using the boundary 
data. □ 

The system $ = (ipk(i)(y), ■ ■ ■ > V'MO^)) °f eigenfunctions and the 
neighborhood {7 in Lemma l£7Tl define a local trivialization of the vector 
bundle V ', namely the map A^y : 7r _1 (£7) — ► [/ x C d , 

-4*,c/(f) = (aj(y))J = i, f G 7r _1 (y), when v = ^a^y^k^y). 

i=i 

Proposition 8.8. The induced bundle structure on dM and the 
Cauchy data set Cq(D) determine a bundle V on M that is bundlemor- 
phic to V . 

Proof: Recall that by Lemma l2~5l we know A r . By Lemma 18771 we 
can find, for any y G M, a neighborhood U and a system $ of the 
eigenfunctions that define a local trivialization in 7r _1 (i7). Thus it 
remains to prove that the boundary data make it possible to evaluate 
the transition functions 

T *,u,*,u :(unu)xc d ^(unu)xc d 
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between any two trivializations A$ t u and A$ 

Let /+ G /- G J 700 , and tel. By Lemma E3 we can find de- 
compositions (|90|) of W+/+(i) and with respect to both Aj,^ 
and fj. Therefore, we can find the representations of the function 

w{y) = Wj + {y,t) + W-f-(y,t), yeUnU 
in the both trivializations A^ ;U and A^^. According to Lemma l8lH 

the function w runs over C°°(M, V), when /+ runs over and /_ 
over J 7 ^. This makes it possible to compare representations 

d d 

w (y) = ^2 a i(v)^k(i)(y) = 5^5j(y)^fc(i)(2/)' y^unu 
1=1 1=1 

when u> runs over C°°(M, V) and, therefore, to determine the transition 
function T^j^- □ 

8.3. Construction of operators on V. We are now in the position 
to complete the proof of the main Theorem 12.91 

Proof: To this far, we have reconstructed (M, g) upto an isometry 
and a bundlemorphic copy V of V. To reconstruct the Hermitian and 
the Clifford module structures and the operators F and D on V, it is 
enough to find their representations in an arbitrary local trivialisation 
A® t u- We will first consider the local trivializations A^ ;U where U C 
M mt and $ consists of eigenfunctions of -Dr,+ only. In this connection 
we will skip index + in the following considerations. 

We start with the reconstruction of the Hermitian structure. Let 
/ G y G U and consider a ball B(y ,r) C U. Denote by /C r = 
JC r (f) the set of those h G J 700 for which 

(92) ((Wf(T)-Wh(T),\ y 5 y )) = 0, for all A, G tT 1 ^), y G B(y ,r). 

Note that condition (j92|) is equivalent to the fact that Wh{T)\ B ^ yo ^ = 
Wf(T)\ B(yo , r) . 

Moreover, it follows from Lemma IQ1 that condition |92l) can be di- 
rectly verified from the boundary data. Therefore, using the boundary 
data we can evaluate 

(93) \\Wf(T)\\l> {B{y0ir)>v) =M \\h\\%. 

As the manifold (M, g) is already reconstructed and W/(T) G C°°(M, V), 
the equation (|93|) makes it possible to find 

(94) \Wf{y ,T)\ yo - hm yol ^ r)) 

for an arbitrary y E U C M mt and / G Again, as W/(T) runs 

over r , ((D r + )°°) when / runs over JF^ , if follows from Lemma l8~H 
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that we can find the Hermitian structure on it 1 (yo) in the trivialization 
A<p u. As yo G U is arbitrary, we can determine the Hermitian structure 
on V | M mt. 

Next we consider projectors Il + : 7r _1 (?/o) — * 7r _1 (2/ ), Ho G M mt . 
Let / G T°° . By Theorem 14.11 and Lemma IQ1 we first find, using the 
boundary data, h G T such that n+W/(T) = Wfc(T). As Wf(y 0l T) 
runs over 7r _1 (?/o) when / runs over JF, this determines Il + in the triv- 
ialization A^jj. Since F = 2U + — I, the chirality operator F is 
determined. 

When U C M int is small enough, there is an orthonormal frame 
(e^y), ... ,e n (y)) on TM\ n . Let y G U andletx(y) = (x 1 ^), . . . , x n (y)), 
x(yo) = 0, be Riemannian normal coordinates centered in y with 
e\ = di at y . By Theorem 14.31 there are generalized sources hji G 
j = 1, . . . , d, I = 1, . . . , n, such that 

Uji(y, T) := Whji(y, T) = x l (y)x(y)(t>k(j)(y), y G U. 

Here \ G C£°(U) is a cut-off function that is equal to 1 in a neighbor- 
hood U C U of y . Moreover, by Lemma IQ1 it is possible to find hji 
from the boundary data. 

Representing the Dirac operator in this local trivialization (cf. fZJ), 

m), 

d d 

Duji(y, T) - A fc (j) x'(y) 0fc(,-)(y) = X] a^?^ li^miy), y eu. 

p=l 

In particular, at y = yo 

(95) Duji(yo,T) = 7(e,)0fc(,-)(j/o)- 

The left-hand side of (j95|) coincides with —iW(d t hji)(y , T). By Lemma 
18.41 the boundary data determines its representation in the local trivi- 
alization A$ [/. Thus we can find a^(yo) G C such that 

7(ej)0fcCj)(Vo) =^aj{(yo)0*(p)(2/o)- 
P =i 

This provides the representation of 7(e/) in the local trivialization A$jj. 
Thus we recover the Clifford module structure in V^int. 

Next, we consider the representation of the Dirac-type operator D. 
For any / G we can find the representations of Wf(y,T) and 
DWf(y,T) = -iW(d t f)(T) in a localjrivialization A$p, U C M int . 
Recall that any functions w G C^°(U, V) can written as w = W/(T) 
with some / G J 700 . Therefore, the set of all pairs 

{{A^,uO / ^f{ T )\u),A^ iU {-iW{d t f){T)\ u )) : / G ^°°} 
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give the graph of the operator D : C£°(U, V) — > C°°(U, V) in this local 
trivialization. This determines D on M mt . 

To recover the Hermitian and the Clifford module structures and 
also the operators F and D near DM, we recall that we have already 
found the bundle structure of V everywhere on M, i.e. the transition 
functions between local trivializations. Consider a local trivialization 
Av# where U HdM ^ and $ = (^k(i)(y), ■ ■ ■ ,^k(d) (y))- By the above 
considerations, we can find the representation of the Hermitian and 
the Clifford module structures and the operators F and D on U D M mt 
in this trivialization. However, these representations are smooth up 
to dM fl U , and we can continue them to DM n U. Thus, we have 
determined the Hermitian and the Clifford module structures and the 
operators F and D everywhere on V. □ 

9. Inverse problem spectral problem. 

9.1 Time harmonic response operator. Let {Af , (fr% |aju}feLi be 

the boundary spectral data corresponding to the self-adjoint Dirac- 
type operator (JHJ), 

Dpu = Du, V(D P ) = {ue H\M, V) : P(u\ dM ) = 0}, 

where P is a (possibly non-local) boundary condition satisfying con- 
ditions (fl9| and (fl8j). Clearly, {4>k\dM}Z=i span the subspace (I — 
P)H 1 ^ 2 (dM, W) and, therefore, determine P. 

We intend to reconstruct, using the boundary spectral data for Dp, 
the fixed-frequency Cauchy data maps, , 



K F X : PL 2 (dM, W) -> L 2 («9M, W), A£/ = u A U- 
is the solution to the boundary value problem, 

Du[ = \u[, P(u{\ dM ) =f, Xia (Dp), 



where a (Dp) is the spectrum of Dp. Af is an integral operator with 
its Schwartz kernel Af(x,y) formally defined as an expansion, 



However, cf. jH^EO] for the scalar Schrodinger operator, this series does 
not converge so that we should regularize it. We use the same method 
as was used in [HO] for the Laplace operator, namely, we represent the 
operator A^ in the form, 

Af/= lim A£/+ f X d v (K p J)dp. 

Observe that in the strong operator convergence in H l l 2 (dM, W), 

9 A Af = hm «9 A Af A , 
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where the Schwartz kernel of d\A^ A is given by 



d x A XA (x,y) = - 



(A-A 7 ) 2 ' 
\\f\<A K 31 

Thus it remains to find rim M _,. +00 A^. 

Lemma 9.1. Let Dp be an operator of form f^j where P satisfies 
conditions ( TT5j) . j TT^) . Then, 

(96) lim Aj/ = / + i 7 (iV)/. 

Proof: Consider the boundary value problem 

(97) Dul = iiiul, P(ul\ dM ) = fePC 00 (dM,W). 

We introduce the boundary geodesic coordinates x(y) = (x'(y),x n (y)) 
where x n = dist(y, dM) and x' are the coordinates of the nearest 
boundary point. Using a special trivialisation of V near dM, see Ap- 
pendix, formulae (I108B and (I109L the Dirac-type operator D takes the 
form 



D = -y(N) (c^ + ^a/,^, D')) , 

n— 1 j 

x n , £>') = - 7 (A r ) ^ 7(e Q ) V Q + ^Z-H(x', x n ) - -y(N)Q(x!, x n ). 

a=X 

Here H(x' ,t) is the mean curvature of the surface dM T = {y G M : 
dist (a;, <9M) = t} at the point y with coordinates (a/, r). 

We look for the solution of |97]) near dM, using a formal WKB 
expansion, 

(98) ul(x',x n )~e-^J2^Mx',Xn). 

Substitution of (|9*£j) to l|9"7}) gives rize to a recurrent system of equations 
with respect to the powers of \x, 

(99) - (i + 7 (A0) u fc + A{x', x n , D')« fc _i + -y(N)d n u k -i = 0, 

where w_i = 0. This system has to be looked at together with the 
boundary conditions, 

' /, for k — 0, 
0, for k jt 0. 



(100) P(u k \ Xn=0 ) 



Equation (|9"9"j) tells us that, for any x, u$ is the eigenspace of 7 (iV) 
corresponding to the eigenvalue — i. We denote the eigenprojectors of 
7 (AQ corresponding to the eigenvalues — i and +i by p~ and p+ 

p; = i(l + vy(JV)), p+=i(l-* 7 (iV)). 



INVERSE PROBLEMS FOR DIRAC OPERATORS 



37 



Applying projector P to (|99|) with k = and using equation j(N)P = 
(I - Ph(iV), we get 

iPu + j(N)(I - P)u = 0. 

Letting x n = in the above equation and using (|100|) , we obtain equa- 
tion flUJ). To justify it, we, however, need to consider further terms in 
the WKB-expansion (l98l) . 

If we apply the projector p~ to (j^Hj) with = 1, we see that 

(101) - id n p~u + p~ A(x', x n , D') p~u = 0, 

where we use that d n r y(N) = and u = p~u . Since 7(e a ) j(N) + 
7(iV)7(e a ) = 0, the differential operator p~ A{x' ,x n , D') p~ is of the 
order 0. Therefore, (|10ip is an ordinary differential equation for p~u 
along the normal geodesic with the initial condition given by (|96|) . This 
determines uq. Next we apply p+ to ([HHl) with k = 1 to obtain ptu\, 

2ip+u 1 = p+ A(x', x n , D')u . 

Together with ()100|) with k — 1, this determines ui\ Xn=0 . Further 
considerations follow the trend above: applying p~ to (|99| with k = 2, 
we obtain an ordinary differential equation along the normal geodesic 
for p~U\\ applying p+ to with k = 2 we determine p+«2 and, by 
(jlOOj) with k = 2, the value W2U„=0i etc. 

Take a finite sum exp(— ftx n ) Xlo^^fc m ijSBl) and cut it by mul- 
tiplying with a smooth function x^n) equal to 1 near dM. Since 
-^(Xlfc=o f i ~ ku k\dAi) = /, and dist (i/x, a{Dp)) > we see that 

z 

IK> - z~^ n ^V~ k UkX\\m{M) < cifT 1 ' 1 . 

implying the desired result on dM. □ 

9.2. From time harmonic data to time domain data. Having Af 
in our disposal, we reconstruct the response operator Ap, see (|35|) . To 
this end we observe that the set 

{u f \ 9 RixR + ■ u f solves ((SSJ) with some / G PC°°(dM x E +) W 7 )}, 

coincides with the Cauchy data set Cq(D). Due to causality, to recon- 
struct Ap it is sufficient to find its values only on PC£°(dM x IR + , W). 

Take / G PC^(dM x R + , W), and let G M+ be such that / = 
for t > . Therefore, u^(t) G V(D P ) for t > T> . Moreover, the graph 
norm of the wave u*(t), 

h f (t)\\l(D) = \\Du f (t)\\h(M,V) + h f (t)\\h {M ,V), 

is independent of t, t > . Since this norm is equivalent to the 
H X (M, I 7 )— norm on V(D P ), we see that 

(102) u f G C , b (M + , H\M,V)), uf\ dMxR+ G C , b {R+, H 1/2 (dM,W)), 
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where C 0> 6(IR + , X) is the class of bounded functions with values in X, 
which are equal to near t = 0. 
Consider the Fourier transform, 



u J 



f(x,X)= / exp(—i\t)u^(x,t)dt. 



By ()102|) . A) is an analytic H 1 (M, V r )-valued function of A in the 
lower half plane, A = v + ifi, \i < 0. For each A, u^(-, A) is the solution 
of the Dirichlet problem (|H7jl with A instead of ifi and /(A) instead of 
/ in the boundary condition. From Af we find u^(X)\au- It follows 
from (fT02l that u f {\)\ qm is bounded along the line {v + ifi : v G M} 
for any < 0. Thus, 



u f \ 9 MxR+ = — exp(-/it) / exp(zz/t)A 1/+ . f(v + ifj,)\ dM dv. 

Summarizing, we see that the induced bundle structure on dM and 
the boundary spectral data {\ k , 4>k\dM}kLi of a Dirac-type operator 
Dp, with P satisfying (fTTl . JTHJ), determine the Cauchy data set, 
Co(D). By Theorem 12.91 this implies that the boundary spectral data 
determine the Riemannian manifold, (M, g), up to an isometry, and the 
vector bundle, V, with its Hermitian and Clifford module structures, 
as well as the chirality operator, F and Dirac-type operator, D, up to 
a Dirac bundlemorphism. Finally, since 4>k\dM, k = 1, 2, . . . , span the 
subspace PC°°(dM, W), we can determine the projector P. Thus we 
can determine Dp. This proves Theorem 12.101 □ 

10. Dirac-type operator in R 3 

In this section we apply our results to the classical example of the 
canonical Dirac-type operator in a bounded domain Mcl 3 . 

Let V = M x C 4 be the trivial bundle over M, each fiber endowed 
with the standard Hermitian inner product of C 4 . The unperturbed 
Euclidean Dirac operator, D has the form 

3 

(103) D u = i ot k d k u + mc 2 a u, u G H 1 (R 3 ; C 4 ). 

k=l 

Here a u , u — 0, ... ,3, are the standard Dirac matrices, 

, k = 1, 2, 3, 

where cx^ are the Pauli 2x2 matrices. Then 

/ 21, for v = 
a^ct M + a^a,, = ^ o for z/ 7^ /i ' ' ^ = ' ' ' ' ' 

and the operator (I103j) falls into the category of the Dirac-type opera- 
tors considered in this paper with F = «oaia2«3 arid 7(efc) = za^. 
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A perturbation Q of a physical nature, consisting of a zero compo- 
nents of the Lorentz scalar and vector potentials, see e.g. jUj, defines 
a perturbed Dirac operator D aq = D + Q. Namely, 

3 

(104) D aq u = i ctk(dk + i(ik(x))u + mc 2 a w + a q(x)u, 

k=l 

where a,k(x),q(x) are C°°-smooth real-valued functions which corre- 
spond to the magnetic vector potential a(x) = a\dx l + ci2<ix 2 + a^dx 3 G 
n x (M) and the scalar potential q(x). Then the potential Q satisfies 
conditions ([TTf and (|12|) and, applying the previous constructions of 
the paper, we obtain the following theorem: 

Theorem 10.1. Let M C M 3 be a bounded domain with a smooth 
boundary. Let D aq and D aq be two Dirac-type operators of form \104\j 
in M . Assume that the Cauchy data sets of these operators are the 
same, i.e. 

C (D a , q ) = C (D a , q ). 

Then q = q and 

h 

(105) a- a = d$ + ^ n m h m , n m E Z, 

m=l 

where {h m } m=1 is a basis of the relative harmonic 1— forms on M, b\ is 
the first relative Betti number of M and <i$ is the exterior differential 
of a function $ G C°°{M,R), $| 9M = 0. 

The 1-forms h m are normalized so that if dM has the external compo- 
nent S and the internal components Ej, j = 1, . . . , b\, then for a path 
7] = 7/([0, 1]) C M, 77(0) G S , and 77(1) G S 3 - we have f v h m = 5 jm . 

Proof: Clearly, the standard Hermitian structure of C 4 , 7(ejfc) = 
and the canonical differential, V 6fe = d^. define the structure of a Dirac 
bundle in M x C 4 with D a q and D a q ~ being Dirac-type operators on 
this bundle. 

By Theorem 12 .9[ there is a bundlemorphism, L : M x C 4 — > M x C 4 
that is a fiberwise isomorphism and, therefore, is given by a smooth 
matrix- valued function L(x) which is invertible for all x G M. Apply- 
ing again Theorem I2.9| we see that D = L~~ X DL, F = L~ X FL, and 
7(e*;) = L~ lr y(ek)L for k = 1,2,3. These imply that a u L = La u , v = 
0, . . . , 3, which means, in turn, that L(x) is a diagonal matrix. As 
the bundlemorphism L preserves the Hermitian structure in each fiber, 
there is a complex valued function k(x) such that 

(106) L(x) — k(x)I, |«(x)| = 1. 

Since the induced bundle structures on the boundary coincide, L(x) = 
I for x G dM. Together with (TTOBll this yields that, cf. |2BJ, k(x) = 
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exp(^(x)), with * G C°°(M,R) such that = on the external 

component of dM and ^(x) = 2nn m , n m 6 Z on the internal m-th 
component of dM. The family of bundlemorphisms fjlOBj) are precisely 
those bundlemorphisms that preserve the bundle structure and, in par- 
ticular, the structure (110411 of the Dirac operator. 

Since D = L^DL, where L is of form (IK)fij) . we see that q = q and 

(107) a=a+dtf. 

We obtain formula (ITOBl) from (HoTI) by writing dW = d$ + ]Cm=i n m h m 
where Q\qm — 0, &i is the first relative Betti number of M and h m is a 
basis of the harmonic 1— forms on M satisfying the relative boundary 
condition $2§§. □ 

Theorem 110.11 is a improvement, for the considered class of pertur- 
bations and domain M C I 3 , of the results in [4~6j . Note that Theorem 
EEI] describes, for the Dirac operator, the relation between the inverse 
problems and the Aharonov-Bohm effect. (For the previous results on 
the Aharonov-Bohm effect for a Schrodinger operator, see e.g. [H3 Elj 
for the Euclidean Laplacian and [2H| for a more general scalar opera- 
tors on a Riemannian manifold.) In physical terms, Aharonov-Bohm 
effect means that not only the the magnetic field da, i.e. curl(a), but 
also the magnetic vector potential a(x) has an effect on the boundary 
measurements. Theorem 110.11 characterizes what part of a besides da 
can be observed from external measurements. 

Appendix. Local Boundary Condition. In this appendix we pro- 
vide a trivialisation of V near the boundary needed in Sections 2 and 
9, and outline the proof that the Dirac-type operator with the local 
boundary condition is self-adjoint. Here, we consider the unperturbed 
Dirac operator 

D T u = D u, V{D T ) = {ue H X (M, V) : P v {u\ dM ) = 0}, 

where D is defined by ((JJ) and the projector Pr by (|26|) . (f28|) . 

First we consider the needed trivialisation of V. In the following, we 
do not change notations for "f(v),v G T X (M), F, etc. when changing 
from one trivialisation to another. 

Let U C M be a coordinate chart near dM with boundary normal 
coordinates, x(y) = (x'(y),x n (y)). We may assume that x{U) — B' x 
[0,e], where B' C W 1 ' 1 is open. Let Wj(x) G 7r -1 (:r), x G Uf) dM, 
be smooth sections such that {iVj(x)}^ =1 is an orthonormal basis of 
7r _1 (x), x G U fl dM. We can continue Wj along the normal geodesies 
so that VatWj = 0, and obtain sections Wj G C°°U, V) that form an 
orthonormal basis for 7r _1 (x), x G U. The map 

^ : ( x, aj(x)wj(x) ) i— > (x, ai(z), . . . , dd(x)) 
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defines a local trivialisation \I/ : 7r _1 (£f) — > U x C d . In this trivialisation 
Vatw is represented by d n u. 

Denoting f) = (x, ip(x, v)), we can define another trivialisation 
f) = (x, u>(x)i[)(x, v)), where in the boundary normal coordinates 

u(x', x n ) = exp f- U 2 1 y #(x', r)rfr^ 

with if (x', r) being the mean curvature of the hypersurface dM T = 
{x G M : dist (x, <9M) = r}. In this trivialisation D has the form 

(108) D = 1 (N){d n + A(x',d x ,;x n )], 

where A(x n ) = A(x', d x >\ x n ) is the first-order differential operator with 
respect to x' which depends on x n as a parameter. Namely, 

n— 1 

(109) A(x / ,9 a/ ;x n ) = - 7 (JV)J] 

a=l 

Here, V a = V 6a . The obtained operator y4.(x', d x r, 0) is called the hy- 
persurface Dirac operator. 

Next we show that Dr fits into the framework of the theory developed 
in [13] implying its self-adjointness and regularity. (For other proofs 
of the self-adjointness see [EH E2|-) To this end, we use the above 
trivialisation $ of V near dM so that D has form (jl08|L fjlOQjt . We 
first show that A = A(r) : L 2 (H) -> L 2 (#) satisfies 

(110) -f(N)A + Aj(N) = and A(0) = A*(0). 

To prove the first identity in (|110|) . we introduce the Riemannian 
normal coordinates on dM T centered at x' Q which are associated with 
the orthonormal frame (ei, . . . , e n _i). Then, 

n—l - 

(111) >y(N)A = - 7 (iV) ^ 7 (iV) 7 (e Q )V Q + ^-#(x', x n ) 7 (iV). 

a=l 

At the point (x' , r), 

n— 1 

V Q (7(A^)u) = -y(N)V a u + -y(y a N)u = -f{N)V a u - ^ S^e^u, 

0=1 

where S a p is the 2nd fundamental form of dM T in the introduced nor- 
mal coordinates centered at x' . Returning to (|111|) . this implies that 

(n—l n—l . \ 

a,P=l a=l J 

By (jSJ, 7(e a )7(e^) +7(e,g)7(e Q ) = — 25 a p and 5 a/ g is symmetric. There- 
fore, the above equation and (llllj) yield the first identity in (jllOjl . 

Similar arguments, employing that in the boundary normal coordi- 
nates, V a ep = S a pN, N = d n , show the second identity in (|110|) . 
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It remains to show that Pr satisfies (fT8j) and (fT9|) . As Pr = h(I + T), 
it follows from (0), © that P r satisfies fifty. To check fify, let P A ps 
be a spectral projector of A(0) satisfying (J2TJ). First equation of 
and © imply that A(I + V) — (I — T) A Thus if A& is an eigenvalue of 
v4(0) with an eigenfunction (f>k, then 

if) k = Pr<Pk - (I - Pr)<Pk 
is an eigenfunction for the eigenvalue — A^. In particular, when A& 7^ 0, 
((<f>k,ipk)) = 0, so that \\Pr<j>k\\ = ~ Pr)<j>k\\ and \\4> k \\ = ||^ fc ||. 
Consider the decomposition, L 2 (dM, W) = Hi © H 2 , where H 2 is the 
0— eigenspace of A(0). The above equations imply that Hi, H 2 are 
invariant subspaces for Pr and 

\\Pf-PapsW = ^= on Hi. 

As H2 is finite-dimensional, then 1 and —1 are not in the essential 
spectrum of Pr ~ Paps- By [IHj, Proposition 3.3, (Pr, Paps) is a 
Fredholm pair, so that Pr is self-adjoint. 
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